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Abstract 



We present a review of the method we have elaborated to compute 
the correlation functions of the XX Z spin-1/2 Heisenberg chain. This 
method is based on the resolution of the quantum inverse scattering 
problem in the algebraic Bethe Ansatz framework, and leads to a mul- 
tiple integral representation of the dynamical correlation functions. We 
describe in particular some recent advances concerning the two-point 
functions: in the finite chain, they can be expressed in terms of a single 
multiple integral. Such a formula provides a direct analytic connection 
between the previously obtained multiple integral representations and 
the form factor expansions for the correlation functions. 
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1 Introduction 



Computing exact and manageable expressions for correlation functions is a central question in 
the theory of quantum integrable models [1-3]. This problem is of great importance, both from 
a theoretical and mathematical view point and for applications to various interesting physical 
situations. Apart from few cases, like free fermions [4-9] or conformal field theories [10], this 
issue is still far from its complete solution. Although several important advances have been 
obtained in the recent years, we are still lacking a general method that could give in particular a 
systematic way to evaluate compact expressions for two-point functions and their long distance 
asymptotic behaviour. 

The aim of the present paper is to give a review of an approach to this problem elaborated in 
[11-13] and in [14-17], together with an account of the more recent progress obtained in [18-20]. 

In our search for a general method to compute correlation functions of quantum integrable 
systems, our strategy is to consider a simple but representative model for which it is possible to 
develop new concepts and tools towards this goal. An archetype of such a model is provided by 
the XX Z spin-^ Heisenberg chain in a magnetic field with Hamiltonian, 

H = H<® - hS z , (1.1) 

where 

M 

# (0) = E KCi + «+i + - 1)} . (i-2) 

m=l 
1 M 

4 = jE°i» [H^,S z }=0. (1.3) 

m=l 

Here A is the anisotropy parameter, h denotes the external classical magnetic field, and cr^ are 
the local spin operators (in the spin-^ representation) associated with each site of the chain. The 
quantum space of states is TL = <8^f =1 7^ m , where H m ~ C 2 is called local quantum space. The 
operators am V ' z act as the corresponding Pauli matrices in the space TL m and as the identity 
operator elsewhere. For simplicity, the length of the chain M is chosen to be even and we 
assume periodic boundary conditions. Since the simultaneous reversal of all spins is equivalent 
to a change of sign of the magnetic field, it is enough to consider the case h > 0. In the 
thermodynamic limit M — > oo and at zero magnetic field, the model exhibits different regimes 
depending on the value of A [1]. The ground state is ferromagnetic for A < — 1, while it has 
magnetisation zero for A > — 1. In the last case the spectrum is gapless for — 1 < A < 1 
(massless regime), while for A > 1 the ground state is twice degenerated with a gap in the 
spectrum (massive regime). 

We are basically interested in the two-point correlation functions of local spins, although 
the results presented here allow us to consider other correlation functions as well. If we restrict 
ourselves to the zero temperature situation, such a problem comes down to the computation 
of the average value of a product of two local spin operators in the ground state \ ip g ) of the 
Hamiltonian 1)1.1(1 : 

g a p(m) = (ip g \a?a^ +1 \i> g ), (a, (3) = (+,-),(-,+), (z, z). (1.4) 

Despite its apparent simplicity, such an object is highly non-trivial to handle. The first problem 
to solve is obviously to determine the ground state \ipg)- A method to diagonalise the Hamilto- 
nian was proposed by Bethe in 1931 [21] and developed later in [22-25]. The algebraic version 
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of the Bethe Ansatz was created in the framework of the Quantum Inverse Scattering Method 
by L.D. Faddeev and his school [26-28]. Different ways to study the correlation functions of this 
model were proposed in the series of works (see e.g. [11, 12, 14-17,29-34]). 

Multiple integral representation for the correlation functions were obtained for the first time 
from the (/-vertex operator approach (also using corner transfer matrix technique) in the massive 
regime A > 1 in 1992 [29] and conjectured in 1996 [30] in the massless regime —1 < A < 1 
(see also [31]). A proof of these results, together with their extension to non-zero magnetic 
field, was obtained in 1999 [11,12] for both regimes using algebraic Bethe Ansatz and the actual 
resolution of the so-called quantum inverse scattering problem [11,13]. In fact, these multiple 
integral representations have been constructed for the elementary building blocks (see SectionEJ), 
since any arbitrary correlation function can be expressed in terms of a linear combination of such 
blocks. One should note however that, although these formulas are quite explicit, the actual 
analytic computation of the corresponding multiple integrals is missing up to now. Moreover, 
the evaluation of the two-point correlation functions (|1.4[1 at lattice distance m is a priori quite 
involved, since the number of terms in the corresponding linear combination of the elementary 
blocks growths exponentially with m (like 2 m ). This makes the problem of asymptotic behaviour 
at large distance extremely difficult to solve in this settings from the present knowledge of the 
elementary blocks. 

In the articles [14, 20] we have derived new multiple integral representations more adapted 
to the two-point correlation functions. One of them [20] is based on the direct re-summation 
of the above linear combinations of the elementary blocks. The second one [14] uses an explicit 
representation for the multiple action of the twisted transfer matrices (see (|2.7|l ) on an arbitrary 
Bethe state. In both cases the number of multiple integrals describing the two-point functions 
(jl.4[) reduces from 2 m to m. 

The development of these methods allowed us to perform further re-summation and to obtain 
representations for the two-point functions on the lattice in terms of a single multiple integral 
[18]. We call this type of representation master equation. The remarkable property of the master 
equation is that it gives a direct analytic link between two general approaches to the computation 
of correlation functions: in the context of the XX Z Heisenberg chain, the first method consists 
in acting with the local operators af and a^ +1 on the ground state ( ip g | to produce a new state 
(ijj(a, /3, m) |, and then in computing the resulting scalar product (ip(a, f3,m) \ ip g ); the second 
method consists in inserting between the two operators of local spin a sum over a complete 
set of states (for example eigenstates of the Hamiltonian) , which gives a decomposition of the 
two-point function in the form 



Using the technique developed in [18, 19], we are able to re-sum completely the form factor 
expansion (|1.5|) and to show that it leads indeed to the master equation obtained by the first 
method. In fact, these two different approaches have a very simple interpretation in the context 
of the master equation. Namely, they correspond to two different ways to evaluate the contour 
integral, by computing the residues in the poles that are either inside or outside the integration 
contour. The first way leads to a representation of the correlation function ( af &^ n -> r i ) m t erms 
of the previously obtained [14] multiple integrals. The second one gives us the form factor type 
expansion of the correlation function (i.e. an expansion in terms of matrix elements of local spin 
operators between the ground state and all excited states). 

This method was generalised in [19] to the time-dependent (dynamical) correlation functions 




(1.5) 



9ap(™,t) = (ip g \ af (0)<t; 



3 

m+l 



{m g ) = {^ g \afe im ai +l e 



-iHt 



(1.6) 
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It is worth mentioning that, up to now, the only known exact results on the dynamical corre- 
lations concern the case of free fermions A = [5,6,9,35-39]. It turns out, however, that the 
methods developed in [14, 18] can be directly applied to the computation of the time-dependent 
correlation functions. In particular, one can obtain a time-dependent analogue of the master 
equation and also a multiple integral representation for g a p{rn, t) in the thermodynamic limit, 
both in massive and massless regime. 

This paper is organised as follows. In Section |21 we briefly recall how to obtain multiple 
integral representations for the elementary building blocks of the correlation functions using the 
algebraic Bethe Ansatz method [11—13], and introduce useful techniques and notations that will 
be used all along the article. In Section |31 we explain how to re-sum these elementary building 
blocks to obtain compact representations for the two-point functions and their generating func- 
tions [14,15]. The problem of asymptotic behaviour for large distances is tackled in Section |1J 
There we consider the toy example of the so-called emptiness formation probability to show 
how the multiple integral representations of Section El can be analysed in the asymptotic limit 
of large distances, both in the massless [16, 17] and massive regimes. We also discuss how the 
methods we have developed could be extended to the case of the two-point functions. Sectional 
and El are devoted to the derivation of the master equation for the correlation functions by the 
two equivalent approaches that have been mentioned above. In the last section we present our 
conclusions and perspectives. 

2 Algebraic Bethe Ansatz and elementary blocks 

Any n-point correlation function of the Heisenberg chain can be reconstructed as a sum of 
elementary building blocks defined in the following way: 

m f 

F m {{e j ,e' j }) = {^ g \J{E e r ei \^ a ). (2.1) 

3=1 

Here | ip g ) is the normalised ground state of the chain and E- 3 ' denotes the elementary operator 

e' e 

acting on the quantum space TCj at site j as the 2x2 matrix of elements E lk ' = 6i !e '6k,e- 

A multiple integral representation for these building blocks was obtained for the first time 
in [29,30]. In this section, we briefly recall how it can be derived in the framework of algebraic 
Bethe Ansatz [11,12]. In general, we have to solve the following successive problems: 

(i) determination of the ground state {ip g \, 

(ii) evaluation of the action the product of local operators on this ground state, 
(Hi) computation of the scalar product of the resulting state with | ip g }, 

(iv) thermodynamic limit. 

The starting point of our method is to use in step (i) the description of the eigenstates 
obtained via algebraic Bethe Ansatz [26,28]. They are constructed in this framework in terms of 
generalised creation and annihilation operators which are themselves highly non-local. Acting 
with local operators on such states in step (ii) is therefore a priori a non-trivial problem. One 
of the key-ingredient of our method, which enables us to compute this action explicitly, is 
the solution of the so-called quantum inverse scattering problem [11, 13]: local operators are 
reconstructed in terms of the generators of the so-called Yang-Baxter algebra, which contains 
in particular these creation/annihilation operators for the eigenstates. Step (ii) can then be 
completed using only the quadratic commutation relations satisfied by these generators [12]. 
The computation of the resulting scalar products in step (Hi) may also present some technical 
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difficulties. In the case of the XXZ Heisenberg chain, it has been solved using the algebraic 
structure of the Yang-Baxter algebra [11,40]. Finally, step (iv) is based on the results of [24,25]. 

Note that this procedure remains essentially the same in the case of the two-point correlation 
functions (see Section |2J). The main difference is that, in step (ii), the reconstruction of the 
corresponding local operators from the solution of the inverse problem gives rise to a more 
complicated combination of the generators of the Yang-Baxter algebra, so that the use of their 
commutation relations to determine their action on the eigenstates involves a more complicated 
combinatoric. 



2.1 General framework 

To compute the elementary blocks (|2.1|) . or more generally any correlation function, the first 
step is to determine the eigenstates of the Hamiltonian and in particular its ground state. 
In the framework of algebraic Bethe Ansatz [26], these eigenstates can be described in terms 
of generalised creation and annihilation operators which are elements of the so-called quantum 
monodromy matrix. In the case of the XXZ chain (|1.1|) the monodromy matrix is a 2 x 2 
matrix, 

with operator-valued entries A,B,C and D which depend on a complex parameter A (spectral 
parameter) and act in the quantum space of states H of the chain. It is defined as the ordered 
product 

T(X) = L M (X)...L 2 (X)L 1 (X), (2.3) 

where L n (X) denotes the quantum L-operator at the site n of the chain: 

r /n_(sinh(H^<) sinhryo-- \ 
^ nlAj ~V sinhr?^ sinh(A - § ' 

Here and in the following, the parameter 77 is related to the anisotropy parameter as A = cosh rj. 

The quantum operators A, B, C and D satisfy a set of quadratic commutation relations 
given by the i?-matrix of the model, and generate the so-called Yang-Baxter algebra. These 
commutation relations imply in particular that the transfer matrices, defined as 

T(A) = trT(A) = A(X) + D(X), (2.5) 

commute for different values of the spectral parameter: [7~(A),7~(/z)] = 0. The Hamiltonian 
(|1.2[) at h = is related to T(X) by the 'trace identity' 

tf (0) = 2 stair/ ^^T _1 (A)| -2Mcoshr/. (2.6) 
(IX 
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Therefore, to diagonalise the Hamiltonian (|1.1[) . it is enough to determine the common eigen- 
states and eigenvalues of these transfer matrices. 

For technical reasons, it is actually convenient to introduce a slightly more general object, 
the twisted transfer matrix 

T K (A) = A(X) + kD(X), (2.7) 
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where k is a complex parameter. The particular case of T K (X) at k = 1 corresponds to the 
usual (untwisted) transfer matrix T(X). It will be also convenient to consider an inhomogeneous 
version of the XX Z chain, for which 

T(A) = L M (A - £ M + v/2) . . . L 2 (X - £ 2 + v /2) L X (A - & + r?/2). (2.8) 

Here, £1, . . . , £m are complex parameters (inhomogeneity parameters) attached to each site of the 
lattice. The homogeneous model IJ1-1JI corresponds to the case where £j = rj/2 for j = 1, . . . , M. 

In the framework of algebraic Bethe Ansatz, an arbitrary quantum state can be obtained 
from the states generated by multiple action of operators -B(A) on the reference state | 0) with 
all spins up (respectively by multiple action of operators C(A) on the dual reference state ( |), 

N N 

|V) =l[B(X J )\0), (if>\ = (Q\l[C(X j ), N = 0,1,...,M. (2.9) 

3=1 3=1 

2.2 Description of the eigenstates 

Let us consider here the subspace 7-i( M / 2 ~ N ) of the space of states TC with a fixed number N of 
spins down. In this subspace, the eigenstates | ip K ({\}) ) (respectively ( Tp K ({\}) |) of the twisted 
transfer matrix T K (n) can be constructed in the form ()2.9|) . where the parameters Ai, . . . , Xn 
satisfy the system of twisted Bethe equations 

^(A i |{A})=0, j = l,...,N. (2.10) 

Here, the function y K is defined as 

N N 

3UHW) = a(ji) 11 sinh(A fc - yu + v) + Kd(n) sinh(A fc - p - 77), (2.11) 
fc=l fc=i 

and a(A), d(X) are the eigenvalues of the operators A(X) and D(X) on the reference state | 0). 
In the normalisation Q2.4|) and for the inhomogeneous model l|2.8|l . we have 

M M 
a(X) = [J S inh(A-^ + r/), d(X) = J| sinh(A - £ a ). (2.12) 

a=l a=l 

The corresponding eigenvalue of T K (fi) on | ^({A}) ) (or on a dual eigenstate) is 

( \<\\\ ( \ TT sinh ( A fc - M + V) , ./ s tt sinh(/i — Afc + ??) 

{A}) = a(/i) M — — tt v— + Kd(ji) \ \ — r-T7 7 \ • (2.13) 

^ smh(A fc - (j,) JLJL sm h(/i - A fc ) 

The solutions of the system of twisted Bethe equations ()2.10|) have been analysed in [41]. In 
general, not all of these solutions correspond to eigenvectors of T K (n). 

Definition 2.1. A solution {A} of the system H2.1U\) is called admissible if 

N 

d(Xj) n sinh(A j - A fc + rj) ± 0, j = 1, . . . , N, (2.14) 

k=l 

and unadmissible otherwise. A solution is called off-diagonal if the corresponding parameters 
Ai, . . . , Atv are pairwise distinct, and diagonal otherwise. 
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One of the main result of [41] is that, for generic parameters k and {£}, the set of the 
eigenstates corresponding to the admissible off-diagonal solutions of the system of twisted Bethe 
equations (|2.1()|) form a basis in the subspace j{( M / 2 - N ) . it has been proven in [19] that this result 
is still valid in the homogeneous case £j = rj/2, j = 1, . . . , N, at least if k is in a punctured vicinity 
of the origin (i.e. < \k\ < kq for kq small enough). Note however that, for specific values of 
k and {£}, the basis of the eigenstates in #M may include some states corresponding to 
unadmissible solutions of (|2.10|) (in particular in the homogeneous limit at K = 1). 

At k = 1, it follows from the trace identity (|2.6|) that the eigenstates of the transfer matrix 
coincide, in the homogeneous limit, with the ones of the Hamiltonian (jl.lj) . The corresponding 
eigenvalues in the case of zero magnetic field can be obtained from (|2.6|) . (|2.L-!|) : 

N 

^ (0) IV'({A}))=^i?(A,)-|V({A})), (2.15) 
i=i 

where the bare one-particle energy E(X) is equal to 

E{X) = sinh(A + ^ s h inh(A-i)- (2 - 16) 
One can similarly define the bare one-particle momentum. It is given by 



2.3 Action of local operators on eigenstates 

e ', e . 

A local operator E?' \ acting in a local quantum space TLj at site j, can also be expressed 
in terms of the entries of the monodromy matrix by solving the quantum inverse scattering 
problem [11, 13]: 

a=l a=l 

This enables us to use the quadratic commutation relations for the generators A,B,C,D of 
the Yang-Baxter algebra to get the action of any product of local operators on arbitrary states 
of the form (JTHJ) [12]: 

N m 

<0|n c ( A *)-n%4(W= E (0|nC(A6), (2-19) 

fc=l j=l 7>C{A} feg-P 

in which the sum is taken over subsets V of cardinal N of the set {Ai,... ,\N+ m }, and the 
coefficients ^p({A}, {e^, e^}) can be computed generically. Thus, the elementary blocks ()2.1|) . 
and more generally any correlation functions (see Section |2J , can be expressed as some sums 
over scalar products of a Bethe state with an arbitrary state of the form ()2.9|) . 



2.4 Scalar products 

We recall here the expressions for the scalar product of an eigenstate of the twisted transfer 
matrix with any arbitrary state of the form (|2.9|) . 
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Let us first define, for arbitrary positive integers n, n' (n < n') and arbitrary sets of variables 
Ai, . . . , A n , m, . . . , [i n and u\, . . . , v n i such that {A} C {v}, the n x n matrix Q K ({A}, {/^[{^j) 
as 



(^«0jfc({A}, = a(p k ) t(\j,[i k ) Y[ sinh(i/ a - fi k + rj) 

a=l 

n' 

- Kd(fj, k )t(/jL k ,Xj) JJ sinh(f a - (j, k - 77), (2.20) 

a=l 

with 

sinhr? 

<( A > M = — T7^ . n / \ r- 2.21 

smh(A — fi) smh(A — /i + rj) 

Proposition 2.1. [11,18,40] Let {Ai,...,Ajv} be a solution of the system of twisted Bethe 
equations \2.W\) . and /Ui, . . . , be generic complex numbers. Then 

N N 

( | [J CQij) | Va{A}) ) = ( ^({A}) | J] B{ H )\ ) 

3=1 3=1 

N N 

n d(X a ) II sinh(/i 6 - A a ) 

= ^ — ^{w tM{X}) J (2 - 22) 

J] sinh(A a - X b ) sinh(/i fe - /u a ) 

iV 

II d(Xa) 

= - ^ detO,({A},M|{A}). (2-23) 

f] sinh(A a - A 6 ) sinh(/x 6 - fi a ) 

a>b 



Remark 2.1. If the sets {A} and {//} are different, the eigenstate | ip K ({X}) ) is orthogonal to the 
dual eigenstate {ip K ({fi}) \. Otherwise 

N 

n d(x a ) 

( MW) I ^({A}) ) = -jj-^ det n K ({X}, {A}|{A}) (2.24) 

n sinh(A a - A&) 

a,b=l 

N 

n d(Xa) f x 

= (-1)"^ det (^(A,-|{A})J . (2.25) 

n sinh(A a - A fe ) 

a,b=l 
a^b 

The equations (|2.22l) — 1|2.25|) are valid for any arbitrary complex parameter k, in particular at 
k= 1. In this case we may omit the subscript k and denote (ip,T,y,£l) = (ipK,T K ,y K ,Q K )\ K=1 . 

Using these expressions for the scalar product and the norm of the Bethe state, one sees 
from equation (|2.19|) that the correlation functions can be expressed as (multiple) sums of 
determinants [12]. 
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2.5 Elementary blocks in the thermodynamic limit 

In the thermodynamic limit, the system of Bethe equations for the ground state turns into a 
single integral equation for the ground state spectral density ptot(A) [25]: 

Ptot(A) + J K{\ - p) p tot (p) dp = ^-t(A, rj/2), (2.26) 
c 

where the contour C, which depends on the value of the magnetic field h, is an interval of the 
real axis in the massless regime and of the imaginary axis in the massive regime (C = [— A^, A/J), 
and the kernel K is given by 

K(X) = isinh27 ? ( 2 . 27 ) 

27rsinh(A + rj) sinh(A — rj) 

For technical convenience, one can also define an inhomogeneous version p(\, £) of this ground 
state density as the solution of the equation 

p(\,0 + I K(\-p)p{p,i)dp = ^t{\i). (2.28) 
c 

Note that /9 to t(A) = p(X,rj/2). In the case of zero magnetic field, this integral equation can be 
solved explicitely and we have 

|A| < 1 : A h = A = oo, (2.29) 
^> = 2( S uJ^- ( y «-*»<>), (2.30) 

A>1: A h = A = -in/2, (2.31) 



More generally, in the limit M — > oo, sums over the solutions Ai, . . . , Ajy of Bethe equations 
for the ground state become integrals over the density p to t- 

N 

^E/W = / Ptot(A)/(A)dA + (l/M), (2.33) 

3=1 c 

for any smooth bounded function /(A). This leads to a multiple integral representation for the 
correlation functions. In particular, the m-point elementary blocks (|2.1|) can be written as a 
m-fold multiple integral of the form [12] 

F m ({e v e' 3 }) = (f[/ dX i) ^({A},{ e „e;})5({A}). (2.34) 

3 

In this expression, the set of integration contours {Cj,j = l,...,m} depends on the regime, 
on the value of the magnetic field, and on the configuration {ej,e'A of the block we consider 
(see [12]). The integrand can be split into two parts: a purely algebraic quantity J^({A}, {ej, e'j}) 
which arises from the commutation relation of the monodromy matrix elements and does not 
depend on the ground state, and a quantity S ({A}) which is the same for all blocks and contains 
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all the informations about the ground state. The latter is actually a functional of the ground 
state density that comes from the thermodynamic limit of the normalised scalar product. In the 
general inhomogeneous case, it is given as 

*«»» = II 8 inhfe- & ) -4f<>^' & "- < 2 - 35) 

l<j<k<m J 

We refer to [12] for an explicit expression of the algebraic part ^({A}, {ej, e' }). Let us just 
mention here that one can essentially distinguish two types of integrals at this level: what we 
will call 'L>-type' integral, that comes from the contribution of the action of an operator D on 
a state of the from (j2.9j) . with its corresponding algebraic part and integration contour Cj = C, 
and the 'j4-type' integral, associated to the action of operator A, with a different algebraic part 
and a contour Cj which is shifted compared to the contour C of the integral equation (j2.2(ij) for 
the ground state density. As the action of operator B is very similar to the successive action of 
operators A and D, it produces in the final result both types of integrals. 

Let us finally note that the representation ()2.34j) for the elementary block (|2.1j) coincides 
exactly for zero magnetic field with the multiple integral representation obtained and conjectured 
in [29,30] (see also [31]) from the q- vertex operator approach, and generalises it to the case of a 
non-zero magnetic field for which the quantum affine symmetry used in [31] is broken. 



3 Re-summation of the elementary blocks 

The method presented in the last section is quite straightforward and gives formally the possi- 
bility to compute any correlation function. However, it has been developed for the computation 
of the expectation values of the monomials T aibl (£i) • • ■ 7a m 6 m (6n)j leading to the evaluation of 
elementary building blocks, whereas the study of the two-point functions involves big sums of 
such monomials. 

Indeed, let us consider for example the correlation function ( c m+1 )• Then, according to 
the solution of the inverse scattering problem (|2.18|) . we need to calculate the expectation value 



m m+1 

< V>({A}) | C(&) • n nta) • B(U+l) ■ II T-\i b ) | V({A}) ). (3.1) 

a=2 6=1 

Since | VK{A}) ) is an eigenstate of the transfer matrix T, the action of YYb^i 1 (S,b) on this 
state merely produces a numerical factor. However, it is much more complicated to evaluate the 
action of Yl™ =2 T{£, a ). Indeed, we have to act first with C(£i) on {ip({\}) | (or with £>(£ m+ i) on 
I VK{A}) )), which gives a sum of states which are no longer eigenstates of the transfer matrix, 
and on which the multiple action of T is not simple. In fact, in the framework of the approach 
of Section the product 0^==2 + ^)(£a) would be computed as a sum of 2 m ~ 1 monomials, 
which eventually would lead to a huge sum of elementary blocks. This is not very convenient, in 
particular at large distance m. Therefore, to obtain manageable expressions for such correlation 
functions, it is of great importance to develop an alternative and compact way to express the 
multiple action of the transfer matrix on arbitrary states or, in other words, to make an effective 
re-summation of the corresponding sum of 2 m_1 terms. 

In this section, we explain two different ways to perform such re-summations. 

3.1 Re-summation with auxiliary integrals 

The results presented in this subsection were first obtained in [14]. We recall here the main 
steps of this re-summation. 
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Let us consider the multiple action of the twisted transfer matrices on an arbitrary dual state 

<oinf=i<?(^), 

N m 
3=1 a=l 

where x\, . . . , x m and fix, . . . , fj,N are generic complex numbers. Using the quadratic commuta- 
tion relations between A, D and C, one can prove 



Proposition 3.1. [14] Let k, xx, ■ ■ ■ ,x m and [ix ■>■■■■> l^N be generic complex numbers. The 
action of Yl™=x %i( x a) on the state ( | Ylj=i ^(h) can ^ e wr itten a s 

N m 

<oin<?o«i)n r «(<> 

3=1 a=l 
min (m,N) 

= e e ^n(K + }>K-}>{^ + },K-})(oi n n c ^)- (3.3) 

n=0 {Ai}={Mc + }U{At a _} «G7+ 6ea_ 

{x}={x 7+ }U{x 7 _} 
I«+I = l7+I=" 

In this expression the set of parameters {/j,} is divided into two subsets {/j,} = {fJ- a+ } U 
such that {/i a+ } n {/Uq,_} = 0. Similarly the set {x} is also divided as {x} = {x 7+ } U {a? 7 _}, 
{x 7+ } n {x 7 _} = 0. These partitions are independent except that #{x 7+ } = = n. 

The sum in V3.$ is taken with respect to all such partitions, and the corresponding coefficient 
R n ({^7+ M x 7- is 9iven by 

R n= \ 1 [ sinh(/x fe - /i ) ] [ sinh(x fe - x a ) ] [ ] [ sinh(/i 6 - /i a ) \ 

^ a>b a<6 a6a + -* 

a,6s«+ a, 667+ 

x Y[ t k( x o,\{x 1+ } U • detO K ({x 7+ },{^ a+ } | {x 7+ } U (3.4) 

aG7- 

The equations (|3.3|) . (|3,4j) are the key-formulae of our re-summation. When applying these 
expressions to particular cases, one obtains directly new multiple integral representations for 
the two-point correlation functions which are essentially different from the ones that result from 
the elementary blocks approach. 

One of the simplest applications of Proposition 13.11 concerns the generating function of the 
two-point correlation function of the third components of spin, which is defined as the expecta- 
tion value 

mK ._ M{A»|Q£ m M{A») 

{Ql ' m) ~ W{A})|V({A})) (3 ' 5) 
of the operator 

m / ., , -, \ l—l 

Qi, 



,m 



n {-^ + -2 1 ■<) =n^)-nm-) • n^ 1 ^, (3.6) 

n=l V J 3=1 j=l j=l 



where | ^({A}) ) is an eigenstate of T(fi) in the subspace 7i^ M / 2 N \ The two-point correlation 
function of the third components of local spins in the eigenstate | ^({A}) ) can be obtained 
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in terms of the second 'lattice derivative' and the second derivative with respect to k of the 
generating function Q3.5JI at n = 1: 



.2 



+ < ^'g^(Q l l,l+m ~ Ql,l+m-l ~ Q?+l,l+m + Ql+X,l+m-l)\ ; (^-7) 



«=1 



Due to the translational invariance of the correlation functions in the homogeneous model, we 
will simply consider the following expectation value: 



(«)in»i«)) 

H 

<v>({a})iv({a}); 



( Q1,m ) = II^feKA}) T^n^m^ • ( 3 - 



In order to evaluate this generating function, one should first compute the multiple action 
of T K (^j) in the r.h.s. of (j3.8j) by means of Proposition 13.11 and then project the result on the 
eigenstate | if)({\}) ) using (|2.22|) for the scalar product. Hereby the expression of the coefficient 

and of the matrices O, Q K can be simplified using Bethe equations for the set {A} and the 
fact that d(^j) = 0. Note also that we can restrict ourselves to the case m < N, since eventually 
we are going to compute the correlation function in the thermodynamic limit. The result can 
be written in the following form 

(Qtm)=f: e n n ^^rr-^^}.^}.^})- ^ 

n=0{A}={A Q+ }U{A Q _}ae 7 - bG 7 + ^ Za) 

{?}={«7+}U{? 7 _} 
I a +I = l7+I=™ 

Here we have combined all the factors in one function T n and extracted explicitly the dependency 
on the subset {£-y_}. We refer to [14] for a more explicit expression. 

Let us now suppose that | ifj({\}) ) is the eigenstate of the inhomogeneous transfer matrix 
which tends, in the homogeneous limit, toward the ground state of the Hamiltonian (jl.lj) . Then, 
in the thermodynamic limit, the sum over the partitions of the set {A} turns, for each given 
n, into an ra-fold multiple integral over the support of the ground state density, just like in 
Section [2.51 for the elementary blocks. As for the sum over the partitions of the set {£}, it can 
be computed in terms of some auxiliary contour integrals. Indeed, it is easy to see that 

e n n tl-n' 1 ' • )) 

{C}=«7-}u{C7 + } a ^- fc e7+ ^ b w 

h+\= n 

n , n m 



J_ I TT dz TT TT sin H z a -£b + v) 
nl J 

r{£} 



n^nn smh^-a) 



n n 

EI sinh(z a - z b ) 
a=l 6=1 

^({z},{X a+ },{X a _}). (3.10) 

11 sinh(z a - Zb + rj) 
a=l 6=1 

Here the contour r{£} surrounds the points £i, . . . , £ m and does not contain any other singu- 
larities of the integrand. Observe that this representation allows one to take the homogeneous 
limit directly by setting £j = rj/2 in the expression. 
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Thus, the sum over partitions in (|3.9j) can be written in terms of multiple integrals. The 
resulting representation for the generating function of the correlation function ( a\ <r^ +1 ) has 
the following form [14]: 



m 1 » n j ~ n m 

«s»>=£sap/nss/*»nn 

n=0 v ; p|^| j=l J c a=l 6=1 



sinh(z a -Zb+r)) sinh(A a - g b 
,2m J I J- sinh(z a - f 6 ) sinh(A a -6 + 77) 

x W n {{X},{z}) -detM K ({X},{z}) •detp(A i ,z fe ), (3.11) 



with 



^({A},{z})=nn 



sinh(z a - A fe + r/) sinh(A b - z a +rj) 
^ ^ sinh(z a - z b + 77) sinh(A a - X b + 77) ' 



and 



(M K ) jk ({\}, {z}) = t(z k , Xj) + k t(Xj,z k ) Yl 



sinh(A a — Xj + 77) sinh(Aj — z a + rj) 



^ sinh(Aj - X a + rj) sinh(z a - Aj + rj) 



(3.12) 



(3.13) 



The integration contour C and the density function p(X,z) are defined in ()2.26fl — (|2.28|) . 

If we had used the expressions of the elementary blocks derived in Sectional we would have 
obtained the generating function (Qf m ) as a sum of 2 m terms, each of them being written 
as a m- multiple integral of the type (|2.34j) . Instead, we have now a representation containing 
only m nontrivial terms. The ra-th term is formulated as a 2n-fold multiple integral, with n 
integrals over the support of the ground state density and n auxiliary contour integrals over 
some auxiliary variables Zj. We will see in Section \5\ that these last integrals play the role of an 
effective re-summation of the form factor series. 

Observe also that, in the homogeneous model, the dependency on the distance m enters 
each integral only as a power of a simple function. This fact might be used for the asymptotic 
analysis of these multiple integrals by the steepest descent method. 

Other two-point functions can be considered in a similar manner. For example, the expec- 
tation value ()3.1|) gives us the correlation function (a^<7~ +1 ). It is clear that one can evaluate 
this correlation function by using first the equations (|3.3|) . (|3.4|) of Proposition I3.1| by acting 
in a second step with the operator B(^ m+ i) on the resulting states, and by finally computing 
the corresponding scalar products via (|2.22j) . All the steps of this derivation are quite similar 
to the ones that we have just described in the case of the generating function ( QI m ). Let us 
merely give here the new multiple integral representation that we obtain by this method for the 
ground-state correlation function ( c m+1 ) in the thermodynamic limit. For simplicity, we 
present the answer in the homogeneous limit and at zero magnetic field: 



m— 1 



1 



n=0 



n!(n + 1)! 



n+l 

n 



dz j 
2ni 



1 



sinh(A n+ i - A n+2 ) 



d n+2 X ( TT sinh ^ a + 2) TT sinh(A a - |) 
V j-j sinh (^ - 1) ' W s inh(A a + 2) 

\ 

■Wn({X},{z}) 



/n+l 

n sinh(A n+ i - z a + rj) sinh(A n+2 - z a ) 

a=l 



n sinh(A n+ i - Xa + ij) sinh(A„ +2 - A a ) 

\a=l J 



x detM K ({X},{z}) ■det [p(X j ,z 1 ),...,p(X j ,z n+1 ),p(X j ,%)] , (3.14) 

7X~ p _L li\\~Zi 



12 



where the contours C and are the same as in (|3.11|) . The analogue Wn({A}, {z}) of the 

function V^dA}, {z}) is 

n n+1 

fl sinh(A a - z b + 77) sinh(z fe - A a + rf) 
W n ({\}, {z}) = — — , (3.15) 

VLJ LJ/ n n n+1 n+1 v ' 

\\ n sinh(A a - h + r]) l[ ]J sinh(z a - z b + rf) 

a=l 6=1 a=l 6=1 

and the (n + 1) x (n+1) matrix M K has the entries 



nC/r \ u \ \ u\ \ TT sinh ( A a ~ A i + v) tt sinh(A i - z b + rj) 

{M K ) 3h = t(z k , A,) - t(X jt z k ) n ginh(Aj _ Xa + 7j) U s[nh(zb _ Aj + ~ } , J<n, (3.16) 

(M«)„+i^ = t(2 fc ,2). (3.17) 
3.2 Alternative method 

There exists another way to reduce the number of terms in the multiple integral representations 
for the two-point functions. In fact, the re-summation which has just been described has been 
performed at the algebraic level: we have computed algebraically the multiple action of the 
twisted transfer matrices on an arbitrary state and, thus, we have avoided any mention of the 
elementary blocks. On the contrary, the method that will be presented below deals directly with 
the elementary blocks in the thermodynamic limit. 

Let us consider again the generating function ( Qf ) for the correlation function of the 
third components of spin. It has been already mentioned that one can, in the multiple integral 
representations (|2.34|) for the elementary blocks, distinguish two types of integrals: the 'D-type' 
integrals (with the original contour C), and the 'vl-type' integrals (with a shifted contour). In 
fact, the generating function (|3.8[) can be decomposed as a sum over elementary blocks obtained 
as expectation values of products of operators A and D only. Such elementary blocks, containing 
only diagonal elementary matrices, can be in general written in the following form: 



F m ({ej, e 3 }) = j d\ x ... j d\ m S({\}) 



c 

x tt sinh(Aj - & + jej - 1)77) sinh(A fc - + (2 - e fc )n) 
sinh(Aj - A fc - (3 - ej - e fc )n) 

where the indexes €j can take two values 1 or 2: ej = 1 corresponds to an '.A- type' integral and 
€j = 2 corresponds to a 'D-type' integral. For simplicity reason, we consider here only the zero 
magnetic field case, but a representation similar to (|3.18|) . with more complicated integration 
contours, can also be written in the case of a non-zero external magnetic field. Hence, the 
generating function ( Q1 m ) can be expressed as a sum of 2 m such terms. It is easy to see that 
the terms which have the same number of 'A-type' integrals exhibit a quite similar structure. 
This observation permits to write the generating function ( Q1 m } as power series on k, 

m 

(Q?, m > = E KS ^( m )' ( 3 - 19 ) 
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where the coefficient G s (m) collects all the terms containing s 'D-type' integrals and m — s 
'-A- type' integrals. This coefficient can be expressed as the following sum, 

G s {m)= Yl F m({ej,ej})- (3-20) 

eiH he™,— m=s 

One can immediately remark from (|3.18|) that the ground state density functional 5({A}) (|2.35f) 
is common for all the terms in this sum. After symmetrisation over the variables A corresponding 
to the integrals of the same type and extraction of the common denominator 



m I 

e^A 1 ,...,A m )=n n sinh(A _ A 



k=lj=s+l J 



n 



sinh(Aj — Afe) 



„, , , , a sinh(A i A/, • /,) sinh(A i - X k - rj) 

x IT sinhCA, - Afc) 

s>j>k>l sinh ( A J - A fc + 7 ?) sinh(Aj - A fc - r?) ' 

we obtain the following representation: 

G s (m)= 1 ! d\ x ... f dA m e^(A 1 ,...,A m )-g s (m,{A}|{e})-5({A i }). (3.22) 

The function Q s {m, {A}|{£}) in l|3.22j) is a rather complicated sum over permutations which 
corresponds to the sum (|3,20|) over all possible configurations of the algebraic part in the expres- 
sion H3.18JI of the elementary blocks. It is possible to express it in a simpler form if we notice 
that it satisfies the four following important properties: 

1. The function G s {m, {A}|{£}) is symmetric under the permutations of the variables £i, 

2. The function e^" 1 " 1 ^' G s {m, {A}|{£}) is a polynomial function of e 2Xj of degree rn — 1. 

3. For m = 1, 

So(l,Ai|£i) =£i(l,Ai|£i) = 1- (3.23) 

4. The function G s ( m , {^}|{£}) satisfies the following recursion relations, 



Hi 



g s (m,{A}|{£})| _ = JJ sinh(A,- - £ a + t?) JJ sinh(A - & + tj) 

aj^k 

x Q s {m- l,Ai,...,Aj_i,A i+ i,...,A m |£i,...,£ fe _i,£ fc+ i,...£ m ), j < s, (3.24) 

m 

g s (m, {A}|{£}), _ = sinh (Aj - £a + r?) [] sinh(A a - £ fc + t?) 

a^k 

x G s -i(m - 1, Ai, . . . , Xj-i, Xj +1 , . . . , A TO |£i, ■ ■ . ,£k-l,£k+l, ■ ■ ■ Cm), j > s. (3.25) 

These properties can be easily proved using the definition of G s ( m , {-M|{£})- They define this 
function in a unique way as for any m they define a polynomial of degree m — 1 in m points. 
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Recursion relations of the same kind as (|H.24|I - (|H.25|) were obtained for the first time by Korepin 
in [42] for the partition function of the six-vertex model with domain wall boundary conditions, 
and the corresponding unique solution was found by Izergin in [43]. The conditions 1-4 are 
very similar to the conditions that characterise the partition function except that they contain 
one more parameter s. However, the expression for the partition function obtained by Izergin 
satisfies these relations for any s. As the solution of the recursion relation is unique, we can 
conclude that the function Q s {m, {A}|{£}) is proportional to the partition function Z m ({A}, {£}) 
and does not depend on s. More precisely, 

Qs{m, {A}|{£}) = — jL- Z m ({\}, {£}), (3.26) 

where the partition function is given by the Izergin formula, 

m m 

n IT sinh(Aj - & + 77) sinh(Aj - 

Zm({A},{£}) = ^S detftA^fc)]. (3.27) 

J] sinh(Aj - A fc ) sinh(£ fc - 
j>k 

We obtain finally the generating function {Q1 m ) as a sum of m + 1 terms, each of them 
being given as a m-fold multiple integral: 

m 

(Q?, m ) = ^« s G s (m), (3.28) 

s=0 

Gs{m) = \ . [ d m X e^(Ax, . . . , A m ) • Z m ({A}|{£}) ■ S({\}). (3.29) 

s\(m — s)\ sinn 77 J 

c 

It is interesting to mention that the first and the last terms in this sum give a representation for 
the emptiness formation probability which will be studied in details in the next section. One 
of the most interesting property of this representation is the presence under the integrals of the 
expression for the partition function of the corresponding six-vertex model with domain wall 
boundary conditions. This is a new and unexpected connection of this very important object 
with the correlation functions of the XX Z spin chain. 

One can note that the two representations (|3,11|) and (|3.28|) of the generating function ( Q1 ) 
that we have obtained in this section are quite different and present different advantages: the 
first terms of (|3.11|) are very simple, but further terms become more and more complicated, 
whereas all the terms of (|3.28j) have more or less the same structure. One can hope that this 
last remark may lead to a common strategy to compute their asymptotics. 



Similar expressions can be obtained for the two-point functions. For example the correlation 

m+l 



function (m) = ( cr m+1 ) can be written as 



ra—l 

E 

s=0 



g + _(m) = 22g + -(m,s), (3.30) 
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in which the coefficients (m,s) are given as the following multiple integrals, 



1 

s!(m — 1 — s)\ sinh m_1 r\ 



~g + -(m, s) „ . , ,„ , ldX 2 ... I dX m I d\+ / cZA_ 



c 



-TT sinh(A. - g fc + rj) sinh(A fc - + v) 
sinh(A_ -Xk + v) 



' m 

JJ 2 smh(A_ - A fc ) 

jij sinh(A + - £ fe ) sinh(A fc - f m+ i) 

sinh(A + - g fc ) sinh(A fc - £ m+1 + rj) \ 
\k=s+2 sinh(A + -A fc -r ? ) ) 

sinh(A+ - Ci) sinh(A_ - ^ + 77) 

■ C3 m _i(A2, • • • , A m J 



sinh( A_ - g fc + ??) sin h(A fc - ^1' 
sinh(A_ — A 

- £k) sinh(A fc - 
sinh(A + - A fc ) 

n 



sinh(A_ — A + ) 

x Z m _!({A 2 , . . . , A m }, {6, • • • , U}) • S({X 2 , ...,X m , X+, A_}). (3.31) 
A very similar representation can be also obtained directly for the two-point function g zz (m). 



4 Towards asymptotic analysis 

We have seen in the last section that it was possible to re-sum, at least partially, the multiple 
integral representation for the two-point function given by the sum over elementary blocks. This 
provides of course a more compact expression but, above all, an expression that seems, due to the 
particular form of the resulting multiple integrals, more suitable for the study of the asymptotic 
behaviour at large distances. In this section, we will see on a simple example how it is indeed 
possible to analyse this kind of integrals. We then discuss the problems that arise when one 
tries to extend this study to either representation (|3.11|) or (|3.28|) of the two-point function. 



4.1 A simple example: the emptiness formation probability 

There exists a particular correlation function for which it is possible to compute the main 
asymptotic behaviour: the so-called emptiness formation probability r(m), which measures the 
probability of formation of some ferromagnetic sub-chain of length m in the (anti-ferromagnetic) 
ground state. It is defined as the expectation value 

r(m) = (V B | n^T^> ( 41 ) 
fc=i 

on the normalised ground state | ip g ) of the chain. Hence, this quantity corresponds to a single 
elementary block, which means that, in the framework of Section [21 it is given as a (single) 
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multiple integral of the type (|2.34j) with m integrations [12,31,44]: 



rim) = lim 1 I 



6,-6n^/2 L J: sinh(£ a - £ 6 ) 



X 



m C j— 1 m 1 

II II sinh(A, - a + V) EI sinh^ -6) 

i=l I k=l k=j+l J 



<r\ j -^^ ^ det[p(A,-,a)]- (4.2) 

c [I sinh(A a - Xb + rj) 

a>b 

Due to its combinatorial simplicity, it has been widely studied recently (see for example [15,45- 
49]). However, the expression Q4.2|) is not convenient for the asymptotic analysis; in particular 
it is not symmetric. Its symmetrised version, obtained in [14], follows directly from the limit 
k — > oo in representations (|3.11j) or (|3.28f) of the generating function ( Q1 m ): 



r(m) = lim — t I d m X TT 



1 



11 sinh(A a - X b + rj) 

a,b=l 

sinh(A a - A 



m 

X 

a<b 



where Z m ({A},{£}) denotes the partition function of the six- vertex model with domain wall 
boundary conditions given by (|3.27|) . From this expression, it is possible to obtain the asymptotic 
behaviour of r(m) using the saddle-point method. This was performed for the first time in [15] 
in the case of free fermions (A = 0), but the method of [15] can be applied to the general case 
as well (see [17] for the study in the massless regime). We briefly recall here the main step of 
this computation and present the result in massless and massive regime. 

To apply the saddle-point method to (|4.3|) . it is convenient to express the integral in the 
following form: 

r(m) = J d m X G m ({X}) e m2 ^(W), (4.4) 
v 

with 



^ m 

S m ({X}) = 2 lo g[ sinh (^ - Xb + rj) sinh(A a - X b - rj)] 

TYI 

a>b 
^ m 

+ — V log[sinh(A a + r,/2) sinh(A a - r}/2)} 



m 

a=l 



1 Hm 1q r/ -2in \™ (detp(Aj,^)) 
■■'-<■ < > LVsinhr]/ J] sinh(£ a - £ b ) 

a^b 



(4.5) 



and 



Gm ( {A) ) Jim j^iM^a . (46) 

det m p(Aj,Cfc) 

In (|4.4|) . the integration domain D is such that the variable of integration Ai, . . . , A m are ordered 
in the interval C = [— A^A^] (i.e. — < Aj < ■ • • < X. m < A^ in the massless case, and 
— zA/j < iXi < • • • < iX m < iAh in the massive case). 
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In the case of free fermions (A = 0), G m ({A}) = 1, and it is easy to see that S m admits 
a unique maximum S m ({\'}) for a set of variables {X[, . . . , X' m } satisfying the system of m 
saddle-point equations: 

dx J S m ({\'}) = 0, l<j<m. (4.7) 

In the limit m — ► oo, the distribution of these variables A's at the saddle point can be described 
by a density function, 

r-W = J™ ^TWY (4 - 8) 



and one can replace sums over the set {A'} by integrals: 



to 

i=i c 



m ^ A' 



/ T^T7P s (A)dA, (4.10) 



, =1 X 'j ~ X 'k m ^°° J A - \ 

3+k C 

for any function / integrable on the contour C. Hence, the system 1)4.70 becomes a single integral 
equation for the density p s (X'), that can be solved explicitely by Fourier transform. Replacing, 
at the leading order in m, the expression of this saddle-point density in the integrals that 
approximate the sums in 1)4.5)1 . one obtains that the main behaviour of the emptiness formation 
probability at the free fermion point in a magnetic field h (\h\ < 4) is given by (see [15] for 
details) 1 

^logr(m) ~ 5 (°) = ilog(i-^). (4.11) 

The general case is slightly more complicated, but follows the same procedure. The main 
problem is that, a priori, we do not know any asymptotic equivalent of the quantity G m (A) 
when m — > oo. Nevertheless, in the case of zero magnetic field, it is still possible to compute the 
asymptotic behaviour of ()4.4)) in the leading order, provided we make the following hypothesis: 
we assume that the integrand of 1)4.4)1 admits a maximum for a certain value A' x , . . . , X' m of 
the integration variables X\, . . . , X m , that, for large to, the distribution of these parameters 
X[, . . . ,X' m can be described by a density function p s (X') of the form 1)4.8)1 on the symmetric 
interval [—A, A] (see 1)2.29)1 . (|2.Hlj) ). and that, at the leading order in m, we can replace the 
sums over the set of parameters {A'} by integrals over this density p s (X') as in ()4.9)) - ()4.1U)) . 

First, like in the free fermion case, it is easy to determine the maximum of the function 
# m ({A}). Indeed, let {A} be solution of the system 

d Xj S m ({X})=0, l<j<m. (4.12) 

In the limit m — > oo, if we suppose again that the parameters Ai, . . . , X m become distributed 
according to a certain density p s (X) and that sums over the Xj become integrals over this density, 

1 For \h\ > 4 the ground state becomes ferromagnetic and the emptiness formation probability is equal to 
(for ft > 4) or to 1 (for h < -4). 
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the system (|4.12j) turns again into a single integral equation for p s , that can be solved explicitely 
in the case of zero magnetic field: 

p s (\) = - y cosfrK) e -2nA 5 ( mass ive case A > 1, ( = -r? > 0), (4.13) 

7T cosh(2nQ 



cosh# 



C\/2cosh^ 



1C 

— — r-, (massless case |A| < 1, ( = irj > 0). (4-14) 



This gives for the maximum of S m ({\}) when m — > oo 2 : 

lim S m ({\}) = - C --f2— S[n ^f} v (A = cosh C >l), (4.15) 

m->oo 2 ^-^ n cosn (2nC) 

n=l v 7 

vr 1 f dj sinh^vr - C) cosh 2 $ „ A . . , 

= log 7 + ^ / . - 2 2 . r^, A = cosC <l . 4.16 

C 2 J uj smh ^ smh ^ cosh 
R-iO 2 2 

The second step is to show that the factor G m ({\}) gives always a negligible contribution 
compared to S m ({\}) at this order in m, at least for any distribution of the variables Xj satisfying 
the previous hypothesis of regularity. Indeed, we can use the integral equation (|2.28|) satisfied 
by the inhomogeneous spectral density for the ground state to express, for any set of variables 
{A}, G m ({\}) in the form: 

Gm{{ }) ~ ftj™.-? det m [ P (A,,&)] ' (417) 

where the kernel K is given by (|2.27|) . If the distribution of {A} is regular enough in the interval 
[—A, A], we can replace, in the limit m — > oo, the integral 

K(\j - fipin^dn (4.18) 
c 

in the determinant by the sum 

1 E^(Ai-A,)^# (4.19) 

where the density function p s (X) describes the distribution of the Xj, j = 1,.. . , m in the limit 
m —* oo. Therefore, 

G m (A) ~ det fe+ ^'^ Y (4.20) 

In the massive regime, this is merely the Fredholm determinant of the integral operator I + K, 
where I denotes the identity operator, and K the integral operator of kernel K l|2.27|l . This 
determinant is given by the infinite product of its eigenvalues: 

oo 

lim G m ({X}) = det(I + K) = 2 TT(1 + q 2n ) 2 , q = e 7 * (massive regime). (4.21) 

n=l 



2 At this main order in m, there exists a unique solution of the integral equation for p s , and we know it 
corresponds to a maximum because S m ({\}) — > — oo on the boundary of 73. 
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In the massless regime, the determinant (|4.'2()|) and its inverse can be bounded via Hadamard 
inequality. Thus, in both regime, we can show that 

lim i-logG m ({A}) =0 (4.22) 

m— >oo ftl 

for any distribution of {A} with good properties of regularity, in particular for the saddle point. 
This means that, at the main order in m, the factor G m ({\}) does not contribute to the value of 
the maximum of the integrand, and that the latter is indeed given by the maximum (|4.15|) - (|4.16|) 
of S m ({\}), p s being identified with the saddle-point density p s . 

Finally we obtain the following result concerning the asymptotic behaviour of r(m) for 
m — > oo (see [17] for the massless case): 

S(°)(A)= lim lQgT l m) , (4.23) 

m— »oo TYl 

= _C_A^ sinhK) (A = cosh C >l), (4.24) 

2 ^ n cosh (2n() v s h K ' 

n=l v ' 

ir 1 f duo sinh % On - Q cosh 2 
= log T + o / . : 1, . ; ,r T^7' (-K A = cosC< 1). (4.25) 



C 2 J u S inh zg- sinh 'S cosh uQ' 

Note that this coincides with the exact known results obtained in [15,47,50] at the free fermion 
point and in [16,46] at A = 1/2, and is in agreement with the expected value in the Ising limit: 

£ (0) (A = 0) = — log 2 (Free fermion case), (4.26) 

S(°)(A = i) = ^log3-31og2, (4.27) 
5 (0) (A) >-oo (Ising case). (4.28) 

A^oo 

Moreover, we can apply the same saddle-point procedure directly at the XXX point A = 1 and 
check that 

5(°) (A = 1) = S (0) (A -> 1+) = S (0) (A -» 1") 

= lQ g( r( r(i)^ ) ra lo s( - 5991 )' ( 4 - 29 ) 

which is in good agreement with the numerical result log(0.598), obtained in [48]. 

In the massless regime, the leading asymptotic behaviour (|4.25|) . was conjectured indepen- 
dently in [51]. In that article was also conjectured the first (power-law) sub-leading correction 
in the form: 

r(m) ~ Am-^e- m2s °, (-K A = cos C < 1), (4.30) 

m^oo 

with 



12 Vvr/ 3(1 - CM 

which is in agreement with the exact results at A = and A = 1/2 (see [47], [16]). It would be 
interesting to check this latter conjecture by analysing corrections to the saddle-point method 
we have presented here. 
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4.2 The two-point functions: attempts and problems 

The long-distance asymptotics of physical correlation functions, such as the two-point functions, 
have attracted long-standing interest. In the case of the XXZ model, some predictions were 
made already a long time ago. 

In the massive regime (A > 1), spin-spin correlation functions are expected to decay expo- 
nentially with the distance and the exact value of the correlation length was proposed in [52]. For 
the XXZ chain in the massless regime (—1 < A < 1), zero temperature is a critical point and 
the correlation length becomes infinite in units of the lattice spacing. The leading long-distance 
effects can be predicted by conformal field theory and the correlation functions are expected to 
decay as a power of the distance. In particular, one expects that, at the leading order, 



A conjecture for the non-universal correlation amplitudes A and A z can be found in [53-55]. 
The exact value of the critical exponents in (|4. 32(1 - 1)4. 331) was proposed for the first time in [56]. 

However, there does not exist at the moment any direct derivation of these predictions from 
the exact expressions of the correlation functions on the lattice. In the last subsection we have 
shown how to determine, at least in the main order, the asymptotic behaviour of the emptiness 
formation probability using the saddle-point method. We could expect to be able to apply the 
same technique to the new multiple integral representation of the two-point function obtained 
in Sectional 

In particular, one can notice immediately that each term of the representation (|3.28|) of the 
generating functional ( Qf m ) has a structure very similar to (|4.3j) . Indeed, it is possible to apply 
to the whole sum a slight modification of the saddle-point technique presented here. It shows 
that, as it should be, there is no contribution of order exp(am 2 ) when m — ► oo. However, to 
obtain the precise asymptotic behaviour of the two-point function, one should be able to analyse 
sub-leading corrections to this saddle-point method, which is technically quite difficult. It is not 
obvious in particular from these expressions that, in the massless regime, the leading asymptotic 
behaviour of the two-point function is only of power-law order. 

If one considers instead (|3.11jl . one can try to make a similar analysis for each term of the 
sum. It can be easily proved that, in the massless regime, the first terms of (|3.11|) decrease as 
powers of the distance. However, it is neither difficult to see that the next terms of the series 
are not negligible with respect to the first ones, which means that one should analyse the whole 
sum to obtain the correct power law asymptotic behaviour. 

5 Complete re-summation for the finite chain 

We have seen in the last section that, although the partial re-summations presented in Section |3] 
contain integrals that can be analysed in the main order via the saddle-point method, the 
asymptotic analysis of the sum itself is much more tricky. It is due to the fact that we have 
to take into account sub-leading corrections to the saddle-point to be able to obtain merely the 
main order asymptotic of either representation (|3.11|) or (|3.28|) . From this point of view, it could 
be more convenient to deal only with one single (multiple) integral instead of a sum, as in the 
case of the emptiness formation probability. 

It is actually possible to re-sum completely representation (|3.9j) to obtain, at the finite chain 
level, what we call the master equation for the two-point function. We will see in particular that 




(4.32) 



(4.33) 
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this master equation sheds a new light on the previous algebraic re-summation of Section 13. II 
by connecting it to the expansion of the two-point function in terms of the form factors of the 
local spin operators. 



5.1 Multiple action of transfer matrices: a complete re-summed formula 

To perform the re-summation of the two-point function, we have used in Section \'S. II the explicit 
formula (|3.3j) for the multiple action of [j^jTsfsa) for an arbitrary set of complex numbers 
{x} on an arbitrary state (0 | Ylf=i [14]. Introducing auxiliary integrals as in 1)3. 10|) . one 

can actually re-sum completely the sum over the partitions of the sets {/i} and {x} in ()3.3|) and 
express this action in the form of a single multiple integral: 

Proposition 5.1. Let k, x\, . . . , x m and . . . , be generic parameters. Then the action of 
lloLi %t( x a) on a state of the form ( | rX/=i C(fij) can ^ e formally written as 

N m f . N , m N 

< o i n cm n t m = t n ^ n ^ {*» • n 

J =1 a=1 r{*}urM j=1 a=1 a=1 



x 



z}) 

N ■ u( _ \ N 

n ZZl-l •y < i«'>-wn'»-<°incw. 

j<k 

where the integration contour T{x} U T{n} surrounds the points 3 x%, . . . , x m and fi±, . . . , //jy an d 
does not contain any other pole of the integrand. 

Proof. The right-hand side in ()5.1|) can be computed by the sum over the residues at the 
points xi, . . . , x m , fj,x, ... , More precisely, it is given by 

min(rfi,iV) 

2 }={x 7+ }u{/i Q _}(I n t e g ranc I), (5.2) 

n=0 {/i}={M a+ }U{/i a _} 
{z}={a; 7+ }U{x 7 _ } 
I Q +I = |7+I=™ 

and one can easily check that the residue of the integrand for {z} = {x 1+ ] U {fJ- a _} is equal to 

^n({z7+M X 7-}'W+}'W-})- D 

To apply this formula to the effective computation of the two-point functions, one should 
give specific values to the parameters x%,..., x m , fJ>i, ■ ■ ■ , /jln- they will correspond in this con- 
text either to admissible solutions of the Bethe equation at k = 1 (parametrising the eigenstate 
in which we compute expectation values) or to inhomogeneity parameters (set to rj/2 in the ho- 
mogeneous case) arising from the reconstruction of local spin operators. To ensure the existence 
of the corresponding contour T{x} U T{n} for these particular values, one has to prove that 
there exists a surrounding of [{xi, . . . , x m } U {fi\, . . . , ^n}] N which does not contain any new 
singularity of the integrand. Since C{z) is a polynomial in z in the normalisation (|2.4j) . such 
singularities could occur only due to the product of the twisted Bethe equations y K (z a \{z}) in 
the denominator. We will see in the following that these poles play a central role in the context 
of master equation. 

3 More precisely, for a set of complex variables {u 1 , . . . ,v{\, the notation T{u} should be understood in the 
following way: is the boundary of a set of poly-disks T> a (r) in C , i.e. r{v} = U„ =1 'D a (r) with ~D a (r) = 

{zeC N : \z k -v a \ =r, k = l,...,N}. 



22 



Solutions of the system of Bethe equations (|2.1()j) were studied in [41] (see also Appendix A 
of [19] for the homogeneous case). Following the arguments of [18], one can easily prove that 
unadmissible and diagonal solutions do not correspond to singularities of the integrand, for 
either detf2 K or (0| Ilj=i C( z j) vanishes in this case. It was moreover proven in [18,19] that, 
for \k\ small enough, all admissible off-diagonal solutions are in the vicinities of the shifted 
inhomogeneity parameters — rj (or of —rj/2 in the homogeneous case) but are separated from 
these points as soon as k ^ 0. For |k| small enough, they are also separated from any admissible 
off-diagonal solution Ai,...,Ajv of the system of untwisted Bethe equations. Thus, just like 
in [18,19], we can formulate the following lemma: 

Lemma 5.1. Let {A} be an admissible off-diagonal solution of the system of untwisted Bethe 
equations, and {£} be a set of inhomogeneity parameters. There exists kq > such that, for 
\k\ < Kq, one can define a closed contour T{\} U r{£} which satisfies the following properties: 

1) it surrounds the points {A} and while all admissible off-diagonal solutions of the 
system \2.1(J\) are outside of this contour; 

2) the only poles which are inside and provide non-vanishing contributions to the integral 
\5.1\) are Zj = A& and Zj = 

3) the only poles which are outside (within a set of periodic strips) and provide non-vanishing 
contributions to the integrand of \5.1)) are the admissible off-diagonal solutions of the system of 
twisted Bethe equations \2.1(l\) . 

Remark 5.1. Lemma 15. II holds also in the homogeneous limit £j = rj/2. 
5.2 Master equation for the two-point function 

Proposition 15. II associated to Lemma 15. II enables us to obtain representations for the two-point 
functions of the finite chain as a single-multiple integral [18]. Let us consider, for a given 
eigenstate | ^({A}) ) of the untwisted transfer matrix, the expectation values 

( a l °"m+l) = ( a l ) + ( a m+l ) ~ 1 + 2 ^72"( ( 5l,m+l - Ql,m ~ Q^m+l + Q2,m)\ J ( 5 - 3 ) 



8k 2 



\k=1 



{of<+i) = lim(<, K < m+ i), = (-,+),(+,-)• (5.4) 

For convenience, we have defined here the following K-deformed spin operators: 

= °? ■ n t ^ ■ n t ^ i (^)' <j = n t ^ ■ ri ■ <>?> 

a=l a=l a=l a=l 

and the operator Q*J k is given by (|3.6|) . Using the solution of the quantum inverse scattering 
problem ()2.18|) and acting with the resulting operators on |, notably by means of ()5.1|) . 

we obtain the following result: 

Theorem 5.1. Let {A} be an admissible off-diagonal solution of the system of untwisted Bethe 
equations, and let us consider the corresponding expectations values i!5.5]) - f5^ in the inhomo- 
geneous finite XXZ chain. Then there exists kq > such that, for \k\ < kq, the following 
representations hold: 



m T t-\ 2ni a-\ T ^™ t\y^a\{z}) 

r{g}ur{A} a_i a_i 
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x , N-l ^ \\T K {^ a \{z}) AT-1 1 

(<7 ^<-+i)= (]v3i)i / n ^ • i+ 2 i • n y K ( Za \{ Z} ) 

r{€}ur{A} 11 r(£ a |{A}) a- 1 

a=l 

^ tt 1 TT 1 det^({A},{z}U{£ m+ i}|{A}) 



X 



L sinh(z j - £ m+1 ) 11 S wh(\j - £i) det^ ft({A}, {A}|{A}) 

E j . "^fo - • II t sinh (A fc - A a + 77) sinh(A fc - A, - 77)] 

a l=i I sinh(A Q -Xp-rj) 



a</3 



a(A /3 ) d(A c 



sinh(A /3 - A a - 7?) fe=1 



Y[ [sinh(A fc - A/3 + 7?) sinh(A fc - A a - rj)] > 
i — i ' 



x det^({4,{A}\{A a ,A^}|{ 2 }) . (5.7) 



1 /" ^dB- U r «(?a|{z}) 7V+1 1 

( ct m ) = T^rrw r II ^i; in II 



(iV + l)! J 11 2*1 ^^alW) 
r{€}ur{A} 11 r(£ a |{A}) a- 1 

a=l 

x tt 1 A 1 »„({*}, {A} U{ft}|{z}) 

|1 sinh(^- - e m+ i) sinh(A j - 6) detjv 0({A}, {A}|{A}) 

(- iV+2 r r s , r s AT+1 

X E { ~ l)a+ " sinYz -Y-ri\ • II [ sinh ^ - *, + V) - % - „)] 



a</3 



a{zp)d{z a ) Y\ r ■ u t ~ _l ^ • ~ M 

| [ [ smh(z fc - zp + rj) smh(zi; - 2 a - 77) J 



sinh(% - 5 a - 77) 11 



x detfi({A},{5}\{5 a ,^}|{A})|, (5.8) 

where we have set (A x , . . . , Ajv+l) = (Ai, . . . , A/y, £1) and (Si,..., ^+2) = (21, • ■ • ,zn+i, Cm+i) 
in J5. 7| ) and m respectively. The integration contours in it5.6)) -i TO)) are sitc/i i/iai i/ie 

on/?/ singularities of the integrand which contribute to the integral are the points £1, . . . ,£ TO and 
Ai . . . , Ajy. 

In the following, this kind of representations will be called master equation. 

Remark 5.2. Expressions (|5.7|) and ()5.8|) are obtained by acting with all the operators on the 
left, i.e. on the dual Bethe state (?/>({ A}) |. It is also possible to act on the right. In that case, 
we obtain for ( af K o"~ m+1 ) a representation that is similar to (|5.7|) . and for {o~\ K o~tm+i) a 
representation similar to (|5.8|) . 

Remark 5.3. Note that the K-deformed two-point correlation functions ( af K o~ K m+1 ), as well 
as the generating functional {Q1 m ), are polynomials in k. They are therefore completely 
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determined by the multiple integral expressions (|5.fij) - (|5.8|) which are valid at least in a vicinity 
of k = 0. The limit k — ► 1 can be reached through analytic continuation. 

Remark 5.4. All these results were formulated in the inhomogeneous case, but hold also in the 
homogeneous limit. We can of course choose | ) to be the ground state of the Hamiltonian 

of the XX Z chain. 

Remark 5.5. From the master equation (|5.6|) . it is easy to come back to the original series 
obtained in Section by decomposing the multiple integral as 

,. N N n N—n 

j n dz 3 = zZ°n j n dz i f n ( 5 - 9 ) 

r{^}ur{A} j=1 n=0 r{£} j=1 r{\} j=1 

Note at this stage that, since the number of poles surrounded by is m and since the 

integrand vanishes as soon as zj = z^, the sum in (J5.9|) is actually restricted to n < m. The 
evaluation of the N — n integrals over the contour T{\} as a sum over the residues leads to a 
sum over the partitions of the set {A} into two subset {A Q _} and {X a+ } of cardinal N — n and 
n respectively. Since the remaining integrals are taken over the contour that surrounds 

only the poles at zj = one can set d(zj) = directly in the integrand. This gives a partial 
resummation, in which we can use the Bethe equation for {A} and take the thermodynamic limit 
as in Section I3TT1 

A similar technique can be applied to (|5.7f) or 1)5.811 . 



5.3 Form factor expansion 

It is also possible, instead of computing the integrals (j5.6)) - ()5.8)) by the sum over the residues at 
the poles inside the contour U T{A}, to evaluate them using the poles outside this contour. 
Due to Lemma 15.1) we know that the only poles that can contribute are the admissible off- 
diagonal solutions of the system of twisted Bethe equations ()2.1U)) . This leads directly to the 
expansion of the correlation functions 1)5.3)1 - 1)5.4)1 over form factors. 

Let us consider for example the two-point correlation function (<r{" c^+i )• Evaluating 1)5.7)1 
by the residues at the poles given by the admissible off-diagonal solutions {u±, . . . , fiN-i} of the 
K-twisted Bethe equations outside the integration contour, one obtains 

m 

n r K (Lm) 

i - + \ _ a=2 

\ a l,K a K,m+l I — 2—1 m+1 

W EI r(Ca\{\}) 

a=l 

y jt i tt i det N n({\},{riu{U+i}\{X}) 

1 = 1 sinh(^ - £ m+1 ) f} x sinh(A i - det^ ^({A}, {A}|{A}) 



X 



E ("^ { ■ - • II -Xa + V) S mh(A fe - A, - n)] 

a b=i < smh (A Q - A/3 - rj) k=1 



a</3 



a{Xj3) d(X c 



FT 1 

Y[ [sinh(A fc - Xp + rj) sinh(A fc - X a - rj)] \ 

fc=l > 



sinh(A/3 — X a — rj) 



x detn K (M,{A}\{A a ,A^}|{//})L (5.10) 
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with (Ai, . . . , Ajv+i) = (Ai, . . . , Ajv, Ci)- I n this expression, we can identify the following matrix 
elements of a± K and cr^ m+l between the Bethe state | ip({\j}i<j<N) ) and the K-twisted Bethe 
state | i> K {{Hj}i<j<N-l) )■ 

m n 

n r(um) n 
0=1 j =i 



( MM) I < m +i I iKW) ) - ^ ^ 

11 t (6j|{A}) n sinli(Aj- - A A ) II sinh(/i fc - Hj) 



a=l j,k=l j,k=l 

j<k j<k 



N-l 

X 



II -T7T^ \ detO({A},Mu{U+i}|{A}), (5.11) 

1A j smh(4 m+ i -/ij) TV 



(^({A})krjv K (M): 



AT-1 

n d(A*i) 
1 j=i 



^(6|M) V ■ W s rT , 

[[ smh(nk - fj,j) [[ smh(A i -A A 

j,k=l j,k=l 
j<k j<k 



N ( N+l 

3=1 K 3 ^ ' K a,/3=l 
a<(3 

N 



f _ a(A a )d(A^) ^ ^ g . nh(Afc _ ^ + ^ g . nh( ^ _ ^ _ ^ 
L smh(A a — Ap — rj) k=1 



a(Xfs) d(X a ) 



N 



Y[ [ sinh(A fc - A/3 + 77) sinh(A fc - X a 



V) 



sinh(A/3 - Aq, - 77) ^ 

x detnCM^A^Aa.A/?}!^})!- (5.12) 

In (|5.12|) . we have defined again (Ai, . . . , Ajv+i) to be equal to (Ai, . . . , A at, Cm+i)- Note that 
the expressions H5.11(l - (|5.12|) of the matrix elements which appear in the summation Q5.1U|) are 
the ones that are obtained when acting on the left with the corresponding local spin operators. 
This is due to the fact that the whole algebraic procedure we have used to derive the master 
equation is based from the very beginning on the principle of action on the left. 
Hence, we obtain 

, _ + . _ V ( V>({A}) I gjg I Mjj) ) • ( MM) I <m+l I ^({A}) ) 

{a ^ a ^+ l} (^({A})|^({A}))-(^(M)|^(M)) ' [ ' ' 

where the sum is taken over all the admissible solutions {fi±, . . . , /J/v-i} of the system of N — 1 
K-twisted Bethe equations. Observe that we did not need to use here the completeness of the 
corresponding K-twisted Bethe states \ ip K ({lJ-})) (see [19,41]) in j{( M /2-N+l) ; as the sum over 
the eigenstates of T K appears automatically as the result of the evaluation of the multiple integral 
(|5.7|) by the residues outside the integration contour. 

Since the resulting expression (|5.10|) - (|5.13ft is again a polynomial in k, we can claim from 
Remark 15.31 that the equality (|5.13|) is valid for all values of k, in particular at k = 1. For 
this value of n and in the homogeneous limit, (|5.13|) represents precisely the form factor type 
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expansion of the correlation function (a± cr^ +1 ), with respect to all the excited states of the 
Hamiltonian. 

Of course, it is now clear that we can proceed in the opposite way, starting from the form 
factor expansion to obtain the master equation H5.7|) and the re-summation of Section [3.11 We 
will discuss this point in more details in the next section. 



6 Dynamical master equation 

The approach described in the previous section can be generalised to the case of time-dependent 
(dynamical) correlation functions. Namely, one can derive a time-dependent analogue of the 
master equations (|5.6|) - (|5,8j) for the two-point functions on the finite chain, 



m+iyL> ' M{a})|^({a}); 



( of (o) ) = xrvL " j/, z; r ■ — , (e.i) 



where | ip{{X}) ) denotes an eigenstate of the Hamiltonian (jl.ip corresponding to an admissible 
off-diagonal solution of the Bethe equations. It is possible to obtain these time-dependent master 
equations either from the formula H5.1|) for the multiple action of the twisted transfer matrices 
or via the form factor approach. Both approaches are described below. 

6.1 Time-dependent generating function 

Let us first consider the dynamical two-point correlation function of the third component of the 
spin in the eigenstate | ip({\}) ). Due to the property [a^, S z ] = 0, we have 

<*i(o)<wi(*)> = mmwm • ( } 

We will explain here how to derive a time-dependent generalisation of the master equation (|5.6|) 
for a generating function of (|6.2[) similar to (|3.5|) . This dynamical generating function is the 
expectation value of the following natural time-dependent generalisation of the operator 4 (|3.6[) : 



Qtmit) = T'- 1 (§) • Tr l+1 (|) • e^ 0) • T~ m (|) • e-^ (0) , (6.3) 
where the k- twisted Hamiltonian is defined similarly as in 

=2sinhr ? ^^T K - 1 (A)| ^-2Mcoshr]. (6.4) 

Using the solution ()2.18|) of the quantum inverse problem and the fact that the transfer matrix 

T K commutes with the twisted Hamiltonian , it is easy to see that, just like in the time- 
independent case, the two-point function (|6,2[) is given as 

( o{ (0) a z m+ i (t)) = 2(<jl(Q))-l + 2D 2 ^{ Ql m (t) ) , (6.5) 



4 In this section we consider the homogeneous chain only, since the local Hamiltonian is well defined only 
in this case. 
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where {Qf m (t) ) is the expectation value of (|6.3j) in the eigenstate | tp({\}) ), 

<0i ~ ( " ) = wumwrn • i6S) 

and D 2 denotes the second lattice derivative defined as in (|3.7|) . 

In order to be able to use Proposition 15,11 to determine the action of the operator H6.3|) on 

the eigenstate | ^({A}) ), one should express the twisted evolution operator e iH * >4 in terms of a 
product of twisted transfer matrices T K . This can be done via the Trotter type formula [57,58], 

^rfWcosM) = & (T K (| + £ ) • T-\T)) ±L , e = i^sinh,. (6.7) 

Using now the identity 

77^/2) = («o(»7/2)d(-r7/2)) "V^-t^), (6.8) 
we arrive at the following representation for the operator Qf m (t): 

QlJt) = Um (2 + e) ■ r r -(|) ■ 7*~(2) . T-(| + e ), (6.9) 

Thus, the problem of evaluating the dynamical correlation function of the third components 
of the spin is reduced to the calculation of the multiple action of twisted transfer matrices on 
the state (^({A}) |. Therefore one can use directly the results of the previous section, with the 
following generalisation of Lemma l5.ll 

Lemma 6.1. Let {A} be an admissible off-diagonal solution of the system of untwisted Bethe 
equations. There exists kq > such that, for < \k\ < hq, one can define a closed contour 
T{A} U T{r]/2} U T{— ry/2} which satisfies the following properties: 

1) it surrounds the points {X}, rj/2 and —rj/2, while all admissible off-diagonal solutions of 
the system h2.1U\) are outside of this contour; 

2) for L large enough, the only poles which are inside and provide non-vanishing contributions 
to the integral \5.1\) are Zj = A&, zj = rj/2 + e and Zj = —rj/2; 

3) for L large enough, the only poles which are outside (within a set of periodic strips) and 
provide non-vanishing contributions to the integrand of \5.1\) are the admissible off-diagonal 
solutions of the system of twisted Bethe equations \2.1U\) . 

Proceeding then to the limit L — > oo, one obtains a time-dependent master equation for the 
generating function (|6.6|) . 

Theorem 6.1. Let {Ai, . . . ,Xn} be an admissible off-diagonal solution of the system \2.1(l\) at 
K = 1. Then there exists kq > such that, for < \n\ < kq, the generating function (Qi m (t) ) 
h6. 6\) in the finite XX Z chain 11. 1\) is given by the multiple contour integral 



1 r N A N 

r{±^}ur{A} j ~ 1 6-1 

JV 

X 



it[E(z b )-E(\ b )]+im{p(z b )-p(\ b )} 



a=i y^\ z a\\ z s) N detAr iZ({A}, {A}|{A}) 
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In this expression, E(\) andp(X) denote respectively the bare one-particle energy and momentum 
\2.1b)) and \2.1T\) ; the integration contour is such that the only singularities of the integrand \5.b\) 
within r{±7j} U T{\} which contribute to the integral are the points {±§1 an d {^}- 

Comparing the dynamical master equation (j6.11j) with the corresponding result (|5.6|) in the 
previous section, we see that the dependency on time appears in the integrand together with 
the bare energy. In the framework of the approach described above, this function arises as a 
Trotter type limit of the eigenvalues of the operators T and T K . The function E(z) has poles 
in the points which explains why the integration contour r{±^} u T{A} in (|6.11l) needs to 
surround not only the point |, but also the point — \. 

Observe also that, unlike in the time- independent case, the equation (|6,11|) holds only in a 
punctured vicinity of k = 0, i.e. for |k| small enough, but non-zero. This is due to two reasons. 
First, the admissible solutions of the twisted Bethe equations are separated from the point — n/2 
only if k 7^ 0. Second, the expectation value {Q1 m {t)) is no longer a polynomial in k: due 

•,rr(0) 

to the presence of the twisted evolution operator e K in (|6.3|) . it has essential singularities 
at K = 0,oo. However, it is easy to see that (Qi m (t)) remains a holomorphic function of 
k everywhere except at these points, which means that the result (|6. llf) can be analytically 
continued from a punctured vicinity of k = to the whole complex plane C*, and in particular 
to the point k = 1. 

6.2 Correlation function ( cr (0) cr + {t) ) 

We have seen in Section E31 that the explicit expansion of the correlation functions over the form 
factors can be obtained from the master equation. It was also mentioned there that one can 
follow an opposite strategy, i.e. sum up the form factors of the local spin operators to derive the 
contour integrals (|5,6|) - H5.8|) , This way to obtain the master equation is more direct and admits 
a very simple generalisation for the time-dependent case. 

Consider, for example, the dynamical correlation function ( c^(0)<7^ +1 (t) ). It would seem 
very natural to compute this correlation function by inserting between the operators crj~ (0) and 
a m+i(t) * ne complete set of the eigenstates of the Hamiltonian (jl.lj) . The point is, however, 
that we have no convenient parametrisation for this complete set. Indeed, it is known that the 
set of the eigenstates corresponding to the admissible off-diagonal solutions of untwisted Bethe 
equations in the homogeneous case is not complete, and that one should take into account 
some unadmissible solutions as well. But, on the other hand, the vectors corresponding to 
unadmissible solutions are ill-defined. 

Nevertheless, we know that the the eigenstates of the twisted transfer matrix T K correspond- 
ing to the admissible off-diagonal solutions of the twisted Bethe equations form a basis of the 
space of states, at least for |k| small enough, but non-zero. We can therefore use these states 
in order to sum up the form factor series. Indeed, let us consider, similarly as in (|5.4|) . the 
K-deformed time-dependent correlation function 

(^(0)< m+1 (t)> = (^e^ 0) -^< m+1 e-^ (0) ^), (6.12) 

with 

<M0)<t+ + i(*)> = l™(%(0)< m+ iW>. (6-13) 

where crj~ K , <7^ m+1 are defined in 1)5. 5 Jl and is given by (|6.4j) . Inserting now the complete 
set of the twisted eigenstates between e ltH * ) ~ lthS ^ an d cr^ m+1 in (|6.12|) . we can immediately 
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write the time-dependent generalisation of the equation (|5.13j) : 

N-l N 

< m°) < m+ iw ) = e- M e n - itE{N) n e ~ uE{Xj) 

{/4 3=1 3=1 

(^({A}) I o7 | MM) > '(MM) I < m +l I ^({A}) ; 



(^({A})|^{A}))-(^(M)|V,(M) 



(6.14) 



where the sum is taken over all the admissible solutions {fix, . . . , fiN-i} of the system of N — 1 
K-twisted Bethe equations. It remains now to repeat all the steps of the corresponding part of 
Section [5] in the opposite order (from (|5,13|) to (|5.7|) ). and we obtain 

iht r N - 1 a„ N ~ 1 N 



< ^(o) < m+1 (t) ) = jf-y f n S n n e 

r{±f }ur{A} j ~ 1 b ~ l 6 - 1 

^1^1 A 1 detiv 0({A}, {2} U { § ) I { A } ) 

X 



a itE(z b )+imp(z b ) TT -itE(\ b )-imp(\ b ) 



X 



n ^(^K4) ' n smh(^ - 2) • n sinh(A . _ 2) detiV n({A}i {A} i { a» 

E C" 1 )^ I ■ "^fo : • n [ sinh ( A * ~ A a + ??) sinh(A fc - A/3 - r?)] 
^ 0= x y. smh(A Q - Xp-rf) k=1 

II [ sinh (Afc - A/3 + 77) sinh(A fc - A a - 77)] j 
k=l ' 



a,/3= 
a</3 



sinh(A/3 - X a - 77) 



x det ^({z},{A}\{A a ,A/3}|{z}) L (6.15) 



iV-l 



In this expression, all the notations are just the same as in (|5.7|) . Note that this result could 
have been obtained also by the method used in Section 16.11 

6.3 Dynamical correlation functions in the thermodynamic limit 

It was explained in Remark 15.51 how one could, starting from the time-independent master 
equations, reproduce the integral representations (|3.11|) - p.l4|) for the two-point functions in the 
thermodynamic limit. A similar method can be applied to the time-dependent case [19], although 
the existence of the essential singularities at ±77/2 in the integrand makes this procedure more 
subtle. 

For simplicity, we consider here only the case of the generating function ( Q1 m (t) ), | ip({X}) ) 
being now the ground state of the Hamiltonian Q1.1JI . Similarly as in Remark 15.51 we can 
decompose the multiple integral in (|6. llf) as 



/N N . n r N 

Y[dZj = J2 C N f H dz 3 f n dz r ( 6 - 16 ) 

r{±2}ur{A} j=1 n=0 r{±f} j=1 r{\} j=n+l 

and evaluate the integrals over the contour r{A}, which leads to a sum over the partitions of the 
set {A} into two disjoint subsets {A} = {A Q+ } U {A a _}. In the time-independent case, we could 
then set d{z) = directly in the remaining part of the integrand, since all the terms proportional 
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to d(z) were holomorphic in the point rj/2. In the time-dependent case, the situation is a priori 
different: due to the essential singularities of the integrand in the points ±rj/2, the contribution 
of the function d(z) in the vicinity of rj/2 (respectively a(z) in the vicinity of —rj/2) does not 
vanish. Nevertheless, since d(z) and a(z) have zeros of order M at z = rj/2 and z = —rj/2 
respectively, one can show that the contributions of the corresponding terms to the total result 
are bounded by C /N\, where C is some constant. Hence, these contributions vanish in the 
thermodynamic limit M,N — * oo, M/N = const. Thus, in the thermodynamic limit, one can 
set d(z) = in the vicinity of z = rj/2 and a(z) = in the vicinity of z = —rj/2. This leads us 
to the following integral representation, in which the integrand is defined differently in the two 
half-planes of the complex plane [19]: 



/n K 1 f J"\ I TT ^fz tt sinh ( Aa ~ Zb ± T. s[nh ( z b -^a + v) 

\ > (n!)2 J j ii 27ri ■ 11 sinh(Aa _ Xb + ) sinh(Za _ Zb + ) 

n=U c T{±2} J a, 0=1 

n 

x Jl ^ E ^- E ^ +im ^- p ^ ■ detM (6.17) 



In this expression, the contour r{±^} surrounds the points ±^ and does not contain any other 
singularities of the integrand. 

Comparing this result with its time- independent analogue l|3.11j) , we see that the determinant 
of densities det[p(Xj, Zk)] has been replaced by a new thermodynamic quantity depending on the 
function 7£^(A, z|{Ai, . . . , A n }, {zi, . . . , z n }). This function is defined differently in the vicinities 
of rj/2 and —rj/2: 



K(X,z\{\},{z}) 



p(X,z), z ~ rj/2; 

sinh(Aft— z+rf) sinh(z— z^+rf) ' 



K~ l n(X z \ n) 11 s ^H^-\+v)smh(z b -z+7 1 ) Z r^-n/2 ^ 6 ' 18 ) 



6=1 



Due to the factors exp(it£ , (z )), the integrand in ()6.17j) has essential singularities in the points 
±7//2. However, in the case t = 0, these essential singularities disappear and the integrals around 
—rj/2 vanish. The remaining part of the integrand has poles of order m at Zj = rj/2. Hence, at 
t = 0, the sum over n in (|6.17|) is actually restricted to n < m, and we reproduce the result of 
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Conclusion 

In this review, we have summarised recent results concerning the computation of correlation 
functions in the XX Z chain by the methods of the inverse scattering problem and the algebraic 
Bethe Ansatz. In conclusion, we would like to discuss some perspectives and problems to be 
solved. 

One of the most interesting open problems is to prove the conformal field theory predictions 
concerning the asymptotic behaviour of the correlation functions. We have seen in Section 2] 
that the new integral representations presented in Section seemed, from the view point of the 
asymptotic analysis, more promising than the original representations in terms of the elementary 
blocks of Section [3 Nevertheless, it is still not clear how to extract the asymptotics directly 
from these multiple integrals. 

A possible way to solve this problem would be to find the thermodynamic limit of the master 
equations. It is natural to expect that, in this limit, one should obtain a representation for the 
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two-point functions in terms of a single functional integral, which could probably be estimated 
for the large time and distance. 

Irrespective of this problem of the asymptotics, the master equation shows a direct analytic 
relation between the multiple integral representations and the form factor expansions for the 
correlation functions. It seems likely that similar representations exist for other model solvable 
by algebraic Bethe Ansatz. It would be in particular very interesting to obtain an analogue of 
this master equation in the case of the field theory models, which could provide an analytic link 
between short distance and long distance expansions of the correlation functions. 

Another interesting further development would be to generalise all these results to the case of 
correlation functions at finite temperature. In this direction, a multiple integral representation 
similar to (|3.11j) was derived recently in [33] for the temperature time-independent correlation 
function, and it would be interested to obtain it in the time-dependent case as well. One 
can also wonder whether there exists a master equation similar to (|5.6|) for the temperature- 
dependent case. It would raise the interesting question of the form factor expansion at non-zero 
temperature. 

It is also well known that, for the case of free fermions A = 0, the dynamical correlation 
functions of the XX Z chain satisfy difference-differential classical exactly solvable equations 
[50,59,60]. It is natural to wonder whether this property holds also for general A, or at least 
for some specific cases. We hope that the multiple integral representations for the dynamical 
correlation functions open a way to study this problem. 
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